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■^j- . interacting. Heat conduction from the hot cluster halo leads to evaporation of the 

relatively cooler galactic halos. The stellar distribution would adiabatically expand as 
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The dynamical evolution of galaxies in clusters is modified if dark matter is self- 



it readjusts to the loss of dark matter, reducing the velocity dispersion and increasing the 
half-light radius. If the dark matter content within that radius was fdm = 25 — 50% of 
the total, as indicated by current observations, the ellipticals in clusters would be offset 
from the fundamental plane relation beyond the observational scatter. The requirement 
that their halos survive for a Hubble time appears to exclude just that range of the dark 
matter cross-section, 0.3 < o~ v jm v < 10 4 cm 2 g , thought to be optimal for reducing 
central halo cusps, unless fdm < 15%. If the cross-section is allowed to vary with the 
relative velocity of dark matter particles, a oc v~ 25 , a new problem of evaporation of 
dark matter arises in the dwarf galaxies with low velocity dispersion. The halos of 
large galaxies in clusters and dwarf galaxies in the Local Group can both survive only 
if 5 < 1.1 or S > 1.8. In either case the problem of central density cusps remains. 



1. Thermal equilibrium in collisional dark matter halos 

Self-interacting dark matter has been proposed recently to alleviate some of the problems 
of galaxy formation arising in the standard cold dark matter models (Spergel & Steinhardt 2000), 
(Wandelt et al. 2000). By allowing for heat conduction into the cooler galactic center, the collisional 
dark matter model reduces the density cusp, which is predicted by numerical simulations in the 
CDM model (Moore et al. 1999), and seems to be in conflict with observations. A drawback of such 
a scenario is that the desired heat conduction takes place in all dark matter-dominated systems, 
producing undesired phenomena in some situations. In clusters of galaxies, the hot cluster halo 
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would tend to evaporate the relatively cooler galactic halos, giving strict limits on the allowed range 
of the self-interaction cross-section. 

Depending on the value of the cross-section per particle mass, a p /m p , collisional dark matter 
can act either as a fluid or as a rarefied gas. The mean free path of dark matter particles is 
A ~ m p /(pa p ), where p is the dark matter density. In the fluid regime, A(r) <C r, the heat diffuses 
out of the hot region (cluster halo) to the cool region (galaxy halo) on the conduction time scale. 
In the opposite optically-thin regime, r < A, the interaction is an elastic scattering of the particles. 

Numerical simulations in both fluid (Moore et al. 2000); (Yoshida et al. 2000) and scattering 
regimes (Burkert 2000); (Kochanek &: White 2000) have shown that in isolated halos the central 
cusp is indeed replaced by a soft core, as desired, for a suitable range of cross-section. However, a 
large range of cross-section appears to be already ruled out by various observations. The experi- 
mental limits and the required rates of core evolution are estimated to narrow the currently allowed 
values to the range 0.1 < cr p /m p < 6 cm 2 g -1 (Wandelt et al. 2000). 

Theoretical arguments seem to be most powerful for the clusters of galaxies. Self-interaction 
tends to make clusters more spherical than follows from the lensing observations. For example, 
(Miralda-Escude 2000) estimated that this constrains the cross-section to be a p /m p < 0.02 cm 2 
g -1 . However, there are serious uncertainties associated with this limit due to statistical and other 
issues. In this paper we consider the effect of thermal evaporation of dark matter halos in clusters. 
Significant evaporation of dark matter halos of the elliptical galaxies would introduce a systematic 
offset in their fundamental plane relations, compared to the field ellipticals. 

In §2 we calculate the dynamics of spherically-symmetric galaxy halos in the fluid regime 
solving the diffusion equation. In §3 we estimate the evaporation time analytically in the scattering 
regime. In §6 we consider the case of a temperature-dependent cross-section that changes the 
dynamics of the inner parts of the halo. Finally, in §5 we show that the evaporation of dark matter 
halos is probably inconsistent with observational scatter of the fundamental plane relation. 



2. Heat conduction in the fluid regime 

In the fluid regime (A < r), the frequent interactions of dark matter particles can be described 
by a diffusion equation 

ds i d 2 or 

at nr z or or 

where S is the entropy, n(r) = p{r)/m p is the number density of dark matter particles of mass m p , 
and T is the temperature defined as 

h B T = X -m p u 2 , (2) 

where u is the one-dimensional thermal velocity of dark matter particles. The coefficient of thermal 
conductivity for the monatomic gas is k ~ c v \nu, where c v = 3/cb/2 is the specific heat per parti- 
cle. From theoretical arguments (Hogan & Dalcanton 2000), the cross-section is either constant or 
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inversely proportional to temperature squared: a = a p (T/T p ) s , 5 = or 2. Also, Wyithe, Turner 
& Spergel (2000) have recently proposed an intermediate case, 5 = 1/2. Therefore, we take 



* m p ' Up V^p/ 

The entropy 5 relates to the equation of state P = nksT oc p 7 with 7 = 5/3. It can be written as 

dS d3 d 1 ! d 3 fc B T 

T ^ = ^ T + P ^n =P A 2^1" (4) 

Assuming that heat conduction is slow compared to the dynamical time in the galaxy, equation (1) 
needs to be solved simultaneously with the hydrostatic equilibrium equation. Appendix A contains 
the detailed implementation of the code. 

As an example, we choose an unspectacular elliptical galaxy NGC 4869 in the Coma cluster, 
situated at a distance of 0.119° from the cluster center (220 kpc for the adopted Hubble constant 
Hq = 65 km s _1 Mpc -1 ). Its effective radius is R e = 4.2 kpc and the half-light luminosity in V 
band is L e = 1.8 x 10 10 L© (Mobasher et al. 1999). The mass-to-light ratio within R e , according 
to the fundamental plane relation for such a galaxy, is 8 in solar units (Gerhard et al. 2001). This 
is about twice the ratio expected for an old stellar population (Binney & Merrifield 1998) (later we 
discuss that the most likely fraction of dark matter mass within R e is approximately 50%), and 
it is consistent with the observed stellar velocity dispersion 200 km s" 1 . We conservatively take 
the total mass of the dark matter halo to be twice the mass contained within the halfdight radius; 
M g = 3 x 10 11 M . 

To set the initial profile for the halo we choose the Hernquist model extending to the tidal 
radius R g> t- 

_ M g R gfi (l + cf 1 
Pg 2ttc 2 r(r + i? 9i0 ) 3 ' {) 

where c = Rg^/Rgfl is the concentration parameter. We consider three values of the core ra- 
dius of the galactic halo, R 9) q = 2.8, 4, 8 kpc. The concentration of the halo may affect its 
(thermo) dynamical evolution. The core radius of about 5 kpc is expected for a galaxy of such 
luminosity from the relation of (Gerhard et al. 2001). 

The density profile of the Coma cluster consistent with X-ray data (Hughes 1989) is 

Pc/ -(l + rV<o) 3/2 ' W 

with po = 10~ 25 hl g cm~ 3 = 2.5 x 10~ 3 M Q pc~ 3 and R c i t o = 290 kpc. The average cluster 
density at the position of NGC 4869 is p d;SCV = 3M d (r)/4vrr 3 = 1.7 x 10~ 3 M pc~ 3 . This tidal 
density determines the truncation radius of the galactic halo, Rg,t> v ^ a Pg,&v ~ 2 Pd,nv an d gives 
Rg tt = 28 (Mg/3 x 10 11 M o ) 1 /3 kpc. 
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The initial profile for the integration is a sum of the halo and cluster distributions. To account 
for the off-center position of the halo in the cluster, we have reduced the central density of the 
cluster profile to the value of the tidal density, po — * Pd,av This is justified as the galaxy size 
is much smaller than the core radius of the cluster and this provides a correct tidal truncation. 
Also, as the galaxy orbits around the cluster center, its tidal radius and other external boundary 
conditions may vary. Using the observed position in the cluster as the average value is reasonable 
since the galaxies are expected to be observed at the location where they spend most of their time. 

Dimensionless equations can be obtained using the characteristic mass M = M g and size 
R = R 9: t of the galaxy halo. The unit of temperature is then To = GMmp/ksR, the unit of 
velocity is vq = (GM/R) 1 / 2 , and the unit of time 



is the optical depth at the tidal radius. 

Figure 1 shows the evolution of the galaxy and cluster halos for the case of constant cross- 
section, 5 = 0. The initial stage, t < 4 to, is driven by internal relaxation in the center but the 
subsequent evolution is accelerated by the heat transfer from the outside. The rate of evolution 
in the cluster is faster than in isolation. As a result of heating, the galactic halo expands and the 
tidal radius shrinks. Eventually, at t = 14.3 to the whole halo evaporates leaving the smooth cluster 
density profile. 

However, the evolution of the temperature profile is not monotonic in time (Figure 2). First, 
the central temperature rises as the heat is conducted inward from the neighboring regions. At 
the same time, the tidal boundary between the relatively hot and cold regions moves inward, the 
halo expands, the confining pressure drops, and as a result the temperature decreases between 1 
and 10 kpc. Then, at t = 12 to the central temperature begins to fall as the whole galaxy becomes 
isothermal. At t = 14 to the central temperature reaches its minimum, even below the initial value. 
At that point the central density is already 100 times lower than initially. Then, the remains of 
the galaxy halo quickly evaporate, and at t = 14.3 to its density falls by another two orders of 
magnitude, while the temperature jumps by two orders of magnitude to match the cluster profile. 

Varying the concentration of the galaxy halo does not change its evolution qualitatively. A 
lower density model (c = 3.5, R g ^ = 8 kpc, same R g ,t) disrupts in 9.9 to, while the higher density 
model (c = 10, R 9: o = 2.8 kpc) takes about 17 to- Thus, the evaporation time is weakly dependent 
on the internal structure of the halo, and the best value can be taken as 14 to- 

The evolution is identical in dimensionless units as long as the galactic halo is in the fluid 
regime. The dark matter cross-section enters only the unit of time, to oc a p /m p . The larger the 



1 

to = Tt tdym 

where t dyn = (GM/R 3 )- 1 / 2 » 1.3 x 10 8 yr is the dynamical time and 



(7) 




(8) 
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cross-section the shorter is the mean free path and the longer is the evaporation timescale, 



(9) 



The stars in the galaxy increase the lifetime of the dark halo by providing an additional 
gravitational attraction in the inner parts. Equation (9) shows that tfi depends essentially on the 
average velocity dispersion of the halo. Adding the fraction of stellar mass, < 50%, would therefore 
increase the evaporation time by about 20%. 

Note, that in the fluid regime the density of the cluster particles does not affect the halo 
evaporation timescale, as long as the cluster represents a much larger heat reservoir. The inward 
propagation of the heat wave is governed primarily by the temperature gradient between the galaxy 
and the cluster. If that gradient is always large, the evaporation time is determined by the internal 
conduction time and is insensitive to the cluster parameters. 



3. Evaporation in the scattering regime 

The situation is different when the particle mean free path is larger than r: in this case the 
fluid equation (1) does not apply. The critical density of the fluid regime, i.e. the average density 
for which the optical depth is unity, is 

»>• = tr <10> 
For G p jm v < 50 cm 2 g _1 , the galactic and cluster halo densities fall below the critical density. If the 
cross-section is considerably smaller, the evolution in the optically-thin regime can be calculated 
at three levels of approximation. 

(i) Cluster particles, which have much higher kinetic energy, would impart a fraction of this 
energy to the relatively cooler galactic particles with each collision. For the same particle mass, 
the energy exchange would suffice to unbind the galactic particle. Thus, each collision with the 
cluster particle removes one from the halo. The mean free path of the cluster particles through the 
galactic halo is A = m p /(a p p g ), and the probability of collision in one crossing time is r/A(r). In 
the optically-thin limit, the collision time per particle at the tidal radius R is then 

\(R) 

*coll ~ Across (R) , (11) 

where the crossing time of the cluster particles is tcross ~ R/v c i ~ (v g/vd) idyri' The collision time 
for all galactic particles is t sc = (p g /pd) t co ii- But since the densities of the galaxy and the cluster 
are similar at the tidal radius, the timescale is of the same order. 

(ii) The calculation can be done more accurately, noting that at any position within the halo 
the rate of evaporation is p = —p(r)/t sc . The amount of bound mass decreases faster because 
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the tidal boundary shrinks as mass is being lost. The resulting mass loss is exponential with the 
timescale roughly t sc /2, depending on the steepness of the density profile. With sufficient accuracy 
we can still take the evaporation time to be t sc . 

(iii) Finally, a heat transfer equation can be used, similar to equation (1), treating the galaxy 
and cluster particles as two separate fluids. In this case the diffusion coefficient has to be modified 
to account for the large mean free path of the particles. The details are described in Appendix B. 

In the scattering regime, the smaller the cross-section the longer is the evolution timescale. At 
the position of NGC 4869 the cluster density is p c \ ps 1.3 x 10~ 3 M pc~ 3 and v c \ ps 10 3 km s _1 . 
This gives for the evaporation time 

Figure 4 shows the evaporation time in the fluid and optically-thin regimes. The fastest 
evolution is achieved for the marginal case A ps R, in which domain we have not made a calculation. 
This intermediate region is sketched in with the dashed line. There, the estimates of the evaporation 
time should agree in the two regimes. For the chosen tidal density, the critical value of the cross- 
section is (J p /m p = 50 cm 2 g _1 . The two time scales, tfi = 4.5 x 10 7 yr and t sc = 7.0 x 10 7 yr, 
agree to within a factor of two. 



4. Other dynamical effects 

Another effect can be important in the fluid regime (pointed out by the referee, Andi Burkert). 
Ram pressure of the cluster halo strips outer parts of the galactic halo as the latter moves around 
the cluster. The restoring gravitational force of the stars balances the ram pressure force on dark 
matter particles only at two half-light radii of NGC 4869. At the tidal radius the gravitational 
force falls to 5% of that. Thus, we expect that on a dynamical timescale the material outside 
about 2 R e ps 8 kpc from the galactic center would be stripped. This would only accelerate heat 
conduction to the inner region and accelerate evaporation of the halo. 

In the optically-thin regime, ram pressure is equivalent to the direct scattering of particles and 
does not add to a new effect. The orbital motion of the halo contributes to the relative velocity 
between the cluster and galactic particles, which is of the same order. Thus, the constraint on the 
cross-section can be strengthened by a factor y/2 (cf eq. [12]). 

Dynamical friction of the galactic orbit in the cluster would also speed up evaporation in either 
regime. As the halo gradually spirals in towards the cluster center, it encounters a higher density of 
cluster particles with higher temperature. However, we ignore this effect as it is unlikely to change 
the constraint on the interaction cross-section significantly. 
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5. Constraints from the fundamental plane relation 

Elliptical galaxies satisfy a tight fundamental plane relation (Djorgovski & Davis 1987). One 
of its projections, the Faber-Jackson relation, has an observed scatter in the central velocity dis- 
persion Alog(j c < 0.07 both in optical bands and in the near-infrared (Mobasher et al. 1999). If 
the ellipticals lose their dark matter halos, their stellar distributions would adiabatically expand 
and reduce the central dispersion, while conserving the total luminosity. This would introduce a 
systematic offset in the Faber-Jackson relation, compared to the field ellipticals that have not lost 
their dark matter halos. 

The adiabatic expansion is characterized by conservation of the azimuthal action, GM(r)r. A 
spherical shell initially at radius r would adiabatically expand to the radius r' such that 

M„(rO r' = [M dm (r) + M„(r)] r. (13) 

The amount of expansion depends on the fraction of dark matter within r : fd m (r) = Md m {r) /M tot (r) . 
The surface brightness of elliptical galaxies is well fitted by the de Vaucouleurs profile. The effective 
radius, R e , characterizes the core of the stellar distribution and contains half the luminous mass, 
M*(R e ) = M*/2. As a result of adiabatic expansion of the core, the effective radius would rise by 
a factor 

^ = I . (14) 

Re 1 - fdm(Re) 

As the core expands, the stellar velocity dispersion falls. If the dispersion profile is flat in 
the core, i.e. o c is determined by the mass distribution within R e , the dispersion would drop by 
Alog<7 c = — log R' e /R e . Then, the galaxy would still satisfy the Faber-Jackson relation within the 
errors if the dark matter fraction is fdm(Re) < 1 - 10~ AlogCTc = 0.15. 

Spectroscopic observations of elliptical galaxies usually have apertures of approximately 3", 
corresponding to 1.5 kpc at the distance of the Coma cluster. The aperture size may be small 
enough compared to R e so that it contains only a small fraction of dark matter, and the measured 
dispersion would respond more weakly to the expansion of the core. As an alternate limiting case, 
we consider that <r c may not change at all, even though the surface brightness does decrease within 
Re- 
in addition to the change of the Faber-Jackson relation, expansion of the stellar core would lead 
to a shift away from the fundamental plane relation. The latter can be expressed as the variation 
with mass of the mass-to-light ratio within R e : M e /L e oc M e - 184 (from Burstein et al. 1997), where 
the mass is defined as M e oc o 2 c R e . The expansion of the core would change the measured value of 
M e while the integrated luminosity L e remains the same. 

Figure 3 demonstrates the effect for NGC 4869. If the core velocity dispersion falls as the 
effective radius expands, o c oc R^ 1 , both the mass and the mass-to-light ratio decrease by the 
factor (1 — fdm) -1 (eq. [14]). If the velocity dispersion does not respond to the core variation at all, 
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both variables increase by the same amount, and the galaxy again moves away from the fundamental 
plane. The magnitude of the offset is determined by the dark matter fraction fdm{Re)- 

The amount of dark matter in elliptical galaxies is harder to measure than in spiral galaxies. 
Fitting results of the X-ray profiles of 42 bright elliptical galaxies (Loewenstein & White 1999) 
require fdm(Re) > 0.2 to explain the spectral parameter j3 = 0.5. The best fitting model with 
Mdm = 3M* has fdm(Re) ~ 0.5. This estimate, that dark matter contributes 50% of the total 
mass of an elliptical galaxy at the half-light radius, is consistent with our model of NGC 4869, 
and also with the model of NGC 4472 obtained using the velocity dispersion of globular clusters 
(from the data of Zepf et al. 2000). We use the value fd m = 0.5 ± 0.25 in Figure 3. 

If the dark matter fraction fdm(Re) > 0.15, elliptical galaxies in clusters need to retain their 
dark matter halos in order to stay on the same fundamental plane as the field ellipticals. The 
fundamental plane relation is observed to be the same in clusters and in the field, and within 
clusters there is no differential radial effect. The evaporation time of the halos, in either fluid or 
scattering regime, must be longer than the Hubble time. Figure 4 shows that this excludes the 
range of the collisional cross-section 0.35 < cr p /m p < 1.1 x 10 4 cm 2 g -1 . 

Note, that similar constraints apply to other types of galaxy. The evaporation of dark matter 
halos of spiral galaxies in clusters would introduce a systematic offset in the Tully-Fisher relation. 
In addition, bulges of spiral galaxies are observed to lie on the same fundamental plane as elliptical 
galaxies. Expansion of the cores of elliptical galaxies in clusters would break this law. 

6. Temperature-dependent cross-section 

The evaporation and heat conduction are modified if the self-interaction cross-section varies 
with the relative velocity of dark matter particles. If self-interaction is coupling- limited 
(Hogan & Dalcanton 2000), the cross-section decreases with temperature as a = a p (T/T p )~ 2 . The 
expected temperature scale for likely candidate particles T p > 1 keV. Other forms of interaction 
are conceivable and therefore we consider a general function for the cross-section varying as a = 
a p (T/T p )~ 5 . Obviously, only the product a p Tp can be constrained. 

Note that different temperatures determine the cross-section depending on the optical depth. 
If r > 1, it is the local temperature of the colder halo; if r < 1, it is the temperature of the hotter 
halo, since it contributes mainly to the relative particle velocity. The optical depth itself depends 
on the temperature, but at the tidal radius the two regimes coincide as T co id — > Thot for r > 1. 

The most optimistic values of the cross-section, as not being in contradiction with observations 
described in the previous section, lie in the optically-thin regime with a p /m p ~ 0.1 cm 2 g^ 1 . Since 
we have the freedom in choosing T p , we fix the temperature scale at u p = (kBT p /m p ) 1 / 2 = 10 3 km 
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s . In this case the evaporation time in the optically-thin regime remains the same: 

/ \ -1 / \ -1+2(5 / \ -1 

t - t ** 1 *(£) (TOT) ( ox^PC-O - < 15 > 

where t>hot = ^c/ and phot = Pd correspond to the hotter halo. The optical depth at the tidal radius 
is also constant regardless of 5; r t = 0.02 a p /m p . Throughout this section, the cross-section a p 
refers to the particular choice of u p . 

In the fluid regime, the diffusion equation is different with 5 > (see eq. [A3]). Compared 
with the constant cross-section case, heat conduction is suppressed in the cooler central parts of 
the galaxy while the heat wave from the cluster propagates faster in the outer parts. Thus, for 
5 > 0, the desired central heating is reduced, but the unwanted evaporation of the outer region is 
increased. 

This requires modification of the external boundary condition in the numerical code. The 
simplest solution is to neglect the scattering events in the cluster halo and move the external 
boundary to the galactic tidal radius. Thus, the galaxy is bound by a fixed cluster pressure at R. 
The evolution proceeds similarly to the case described in §2, but somewhat faster in the outer parts. 
The evaporation time is tfi = q(5)to(5), where q(5) is decreasing with increasing 5; q(5) w 1 — 10. 
However, the unit of time to is scaled up with 5. Thus, the evaporation time in the fluid regime is 

t ^3.0x 10 ^)(^)( I ^)'- 2 V, da) 

where f co id = {GM / R) 1 / 2 is the unit of velocity of the colder halo. For 5 = 2 and f co id = 215 km 
s _1 , appropriate for NGC 4869, tfi ~ 3 x 10 7 (cr p /m p ) yr. The fluid regime applies to NGC 4869 
for a p /m p > 50 cm 2 g _1 . 

Even if we take o p jm p = 0.1 cm 2 g -1 , so that the optical depth is very low at the tidal radius 
of NGC 4869, the central part may become optically-thick if 5 > 0. Because of the lower velocity 
dispersion, the cross-section would rise steeply there. For example, for 5 = 2, the inner 8 kpc arc 
optically-thick. Because of the negative heat capacity of the self-gravitating central region, heat 
transfer into the center would only lower the temperature and increase the cross-section. A natural 
lower limit of the temperature is the velocity dispersion provided by the stars, of order 200 km s _1 . 
The central density of dark matter would fall until the optical depth falls below unity. If we require 
that the dark matter mass within r s does not fall by more than 50%, the cross-section must not 
decrease with temperature faster than 5 « 1.1. 

Note that for 5 > the optical depth increases for smaller halos, by a factor (u p /vh_ ot ) 2S (cf 
eq. [8]). In the optically-thin regime, the cross-section increases by the same factor. Therefore, 
the halos of dwarf galaxies would evolve faster within the Local Group (t>hot ~ 220/\/2 « 150 km 
s _1 ) than within large clusters («hot ~ 1000 km s _1 ). For example, Carina, Draco, and Ursa Minor 
are entirely dark-matter dominated dwarf galaxies, with the central velocity dispersion w co id ~ 6.6 



February 1, 2008 



10 



km s _1 (Mateo 1998). If the dark matter particles in their halos become optically-thick to self- 
interaction the cross-section grows even further, by a factor {u p /v co \&) 2& compared to the (5 = 
case. For the large or small enough values of 5, depending on the regime, the dwarf halos would 
completely evaporate. This allows us to constrain a possible range of 5 from the requirement that 
the dwarf halos survive for the Hubble time. 

The density of the Galaxy halo at a distance of 70-100 kpc, where the dwarf spheroidals are 
located, is of order p^ot ~ lO -4 M pc -3 . The estimated masses of the dwarf halos are M « 
2 x 1O 7 M , and the corresponding tidal radii R « 2 kpc. We shall fix the normalization of the 
cross-section at <J p /m p = 0.1 cm 2 g~ x , suitable for large galaxies. Then, the dwarf halos would 
satisfy the condition n < 1 as long as 5 < 2.8. In fact, in the cores of these galaxies the density is 
much higher than average, po ~ O.4M0pc -3 within the core radii ro ~ 200 pc (Mateo 1998). In 
order for the centers to remain in the optically-thin regime, the halos need 6 < 1.7. From equation 
(15), the evaporation time is longer than 10 10 yr if S < 1.5. Thus, in the scattering regime 5 should 
be small enough to rescue the dwarfs. 1 

If the slope 5 is high enough (5 > 1.7), the dwarf galaxies may fall into the optically-thick 
regime even for u p /m p = 0.1 cm 2 g _1 . From equation (16), the evaporation time is longer than 
10 10 yr if 5 > 1.8 (or 5 > 1.5 for o p /m p = 1 cm 2 g _1 ). In other words, the cross-section should 
increase substantially enough that heat conduction slows down. 

Thus, the presence of the dark matter-dominated dwarf galaxies in the Local Group rules out 
a certain range of the temperature slope, 1.5 < 5 < 1.8, even for the most optimistic value of 
a p /m p = 0.1 cm 2 g" 1 . Together with the constraint from NGC 4869, this excludes the values 
l.K 5 < 1.8. 

7. Conclusions 

We have considered astrophysical constraints on the possible cross-section of self-interacting 
dark matter. The fundamental plane relation of elliptical galaxies requires that their dark halos 
largely survive the heating from the hot cluster halos. This condition applies to galaxies with the 
half-light dark matter fraction > 0.15. For a velocity-independent cross-section, the evapora- 
tion time of such halos is shorter than the Hubble time if 0.3 < cr p /m p < 10 4 cm 2 g _1 . 

The parameter space can be enlarged by assuming a p /m p = 0.1 cm 2 g -1 and allowing for 
a velocity dependence of the cross-section. The large halos would be unaffected, as long as they 
remain in the optically-thin regime, but the dwarf halos within larger halos would evaporate if the 
optical depth rises towards unity. The presence of dwarf spheroidals in the Local Group excludes 
the range of the power-law slope 1.1 < 5 < 1.8. 



1 These constraints become more severe if we take a larger normalization of the cross-section, a p /m p = 1 cm 2 g 
In this case, the cores of the dwarfs remain optically-thin for 8 < 1.1 and survive for the Hubble time if 8 < 0.9. 
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Combined with other constraints in the literature (Wandelt et al. 2000), this appears to ex- 
clude the most interesting range of the cross-section for the purpose of lowering the core densities 
of galactic halos. In order to reduce the central density cusp by scattering, the optical depth 
at the scale radius of the halo R g fl should be close to unity (Spergel & Steinhardt 2000). For a 
Hernquist model of the halo, r{R g fi) exceeds the optical depth at the tidal radius, r t , by a factor 
(1 + c) 2 /4 ~ 10 — 30. We can therefore expect the optimum regime to be Tt ~ 0.03 — 0.1, but as 
Figure 4 demonstrates, this regime would lead to evaporation of the halo in less than a Hubble 
time. 

Thus, the original motivation for considering self-interacting dark matter may require colli- 
sional cross-sections too large to be consistent with observations of galaxies in clusters. However, 
the cross-section in the optically-thin regime {<J p /m p ~ 10 _1 cm 2 g _1 ) may still have very interest- 
ing consequences for the growth of central black holes (Ostriker 2000). 

We would like to thank Martin Rees, David Spergel, and Paul Steinhardt for valuable discus- 
sions, Nick Gnedin for advice with the numerical scheme, and the referee Andi Burkert for useful 
comments. This work was supported in part by PPARC and by NSF grants ASC-9740300 and 
AST-9803137. 



A. Numerical scheme for the fluid regime 

For computational convenience, the equations are rewritten in Lagrangian form using the new 
variable of integration dM r = A.-Kpr 2 dr. In dimensionless form they become 



dr 1 



8M r 4vrpr 2 ' 
dP M r 



8M r 4vrr 4 ' 



(Al) 
(A2) 



9 dtp^'dMrX 9 dMr)' [ ' 

The unit of time (eq. [7]) absorbs all parameters of the problem. The equation of state is enforced 
by defining the specific entropy as P = sp 7 . The central differencing scheme is 

rf+r - rt a ' 3-^1/1 1 



mi+i-rrii An 2 yp^ i Pi+i^+i 

-l-rrij l fruj m 
4vr 2 I rf r 



= 0, (A4) 

7 7 m i+1 -mil f rrii m i+1 \ 
SipJ - s i+1 pj +1 — -4 +-4— =0, (A5) 



+1 



T i z2± = r = Fi+1/2 - Fi ~ 1/2 (AG) 

At 1 (m m -m;_i)/2' 1 ; 
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where T* is the solution at the previous time step, and 

j? 1 / 4 ^5+1/2 . 4 rpS+l/2\ T i+ i - Ti 

Fi+i/2 = ^{r i piT i + r i+lPi+1 T i+1 J — — — , (A7) 

_ Inpj+i - In pi 

tti = : . A8 

In r i+ i - In n 

The time- variable local slopes of the density distribution, «j, provide the most accurate integration. 
The grid is logarithmic in radius. The inner boundary condition assures constant mass of the 
first mesh point, p\r\ = const. The outer boundary condition fixes the external pressure of the 
continuous cluster distribution, pnTn = const. 

The finite difference equations are solved using the iterative correction procedure. Corrections 
to the previous configuration are sought: 5r = r — r*, bp = p — p* , 5T = T — T*. First, the 
temperature is updated using equation (A6). Then the entropy is fixed, 

s,s w?W < A9 > 

and 5 P i are calculated from i = N downward and 5ri are calculated from i = 1 upward. The cor- 
rections are added to the previous solution and equations (A4,A5) are iterated until the corrections 
converge. The time step is chosen such that the relative corrections 5r/r, 5p/p, ST/T do not exceed 
10%. 



B. Heat transfer in the optically-thin regime 

The coefficient of thermal conductivity for the monatomic gas is n ps c„nA 2 /t co ii, where c v = 
3kB/2 is the specific heat per particle. In the optically-thick regime, A ps ut co \\ and k = c v \nu. In 
the optically-thin regime, A(r) > r, this should be modified to replace A, which is now larger then 
the size of the system, with titdym the characteristic distance traveled by the particles (similar to 
heat conductivity in globular clusters; (Lynden-Bell & Eggleton 1980)). Thus, the conductivity is 
multiplied by the local optical depth r(r) = t& yIl (r) /t co \\{r) squared. Taking t^yn = (47rGp) -1 / 2 
and t co n = nip/ {pa p u) , we have 

_3 *#V/g fg p V k B T P 
K °-2 m ^ ap X {mJ 4nGm p (m) 

For simplicity, we assume the density and temperature of the cluster particles to be fixed and 
uniform throughout the halo. The change of the entropy of the galactic halo particles is the sum of 
the internal heat conduction with other halo particles and the direct heating by cluster particles. 



dS dQ cond dQ heat 1 d 2 dT T t 



at at at nr z or or i, r 
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With a new unit of time, to = 3/(647r 3 r() t^ yD , the equation can be written in dimensionless form 
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where t' = 16ir 2 



/(PclVd)- 
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Fig. 1. — Evolution of the density profile of NGC 4869 in the fluid regime. The galaxy halo is 
initially a Hernquist model with the concentration c = 7, the cluster is similar to an isothermal 
sphere. The origin of the coordinate system is placed at the galactic center which does not coincide 
with the cluster center. Dashed vertical line marks the initial tidal radius of the galaxy halo. The 
outputs are at the times indicated, in units of to = 6.5 x 10 4 (a p /m p ) yr. As time progresses the 
galactic density falls continuously and approaches the ambient cluster profile. 
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Fig. 2. — Evolution of the temperature profile of NGC 4869 in the fluid regime. The temperature 
is converted to the velocity units, independent of the particle mass m p . 
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Fig. 3. — Fundamental plane relation for the elliptical galaxies (solid line) and its scatter (dots). 
The half-light mass-to-light ratio M e /L e is in solar units. If NGC 4869 loses its dark matter halo, 
it would deviate from the plane either upward (if the core velocity dispersion stays constant as the 
effective radius expands) or downward (if a c oc -R" 1 ). The length of solid arrows is calculated using 
fdm = 0.25, a minimum dark matter fraction within R e , while dashed arrows are for best-fitting 
model, f dm = 0.5. 
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Fig. 4. — Evaporation time of the dark matter halos of elliptical galaxies in clusters in the fluid 
(t< > 1) and scattering (tj < 1) regimes. The lines broke into dashes in the marginal range, for the 
optical depth at the tidal radius 1/3 < r t < 3. The arrows indicate the range of the cross-section 
excluded by the requirement that the halos survive for 10 10 yr. 



